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In  a  split  beam  sonar  system  the  signal  to  noise  ratio  an-1  the 
azimuthal  resolution  may  be  altered  by  varying  the  configuration  of  the 
subarrays.  This  paper  analyses  the  performance  of  such  a  sonar,  taking 
account  of  the  noise  correlation  between  the  subarrays  to  optimise  the 
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THE  OPTIMISATION  OF  THE  ARRAY  CONFIGURATION  FOR 


A  SPLIT  BEAM  SONAR  SYSTEM 


I.  INTRODUCTION 

A  common  requirement;  in  the  design  of  a  sonar  is  to  achieve 
a  high  azimuthal  resolucion  subject  to  constraints  on  the  sonar  array 
size.  One  method  of  doing  this  is  to  use  a  split  beam  technique  in 
which  the  array  is  split  into  two  subarrays  which  both  ensonify  the 
sonar  target.  The  target  resolution  is  then  determined  by  measuring 
the  time  difference  between  the  signal  returns  for  the  two  arrays. 

The  measurement  of  the  time  difference  is  usually  carried  out  by 
measuring  the  phase  difference  between  the  signals  received  at  each 
subarray. 


For  a  given  overall  array  size  the  performance  of  such  a 
sonar  may  be  varied  by  altering  the  dimensions  of  the  two  subarrays. 

When  this  is  done  the  noise  measured  at  each  subarray  may  not  generally 
be  assumed  independent.  This  paper  analyses  the  performance  of  a  split 
beam  sonar  under  these  conditions  to  determine  the  optimum  configuration 
of  the  subarrays. 
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2 .  THE  AMPLITUDE  AND  PHASE  AT  THE  SUBARRAYS 

Suppose  that  a  split  beam  sonar  system  consists  of  two  sub¬ 
arrays  of  width  b  and  a  spacing  between  array  centres  of  d.  The  signal 
received  at  the  arrays  is  resolved  to  a  number  of  range  cells.  Then 
for  each  range  cell  the  amplitude  of  the  received  signal  is  calculated 
by  adding  the  amplitude  received  at  each  sub-array,  and  the  phase  is 
calculated  by  subtracting  the  phase  received  at  each  subarray.  It  will 
be  assumed  that  the  phase  calculated  for  each  subarray  is  referenced 
to  the  reflected  target  signal  received  at  the  electrical  centre  of 
the  total  receiving  array.  This  point  is  marked  0  in  Fig.  1.  It  is 
assumed  that  in  any  range  cell  the  signal  reflected  from  the  target  may 
be  represented  by  a  reflection  from  one  azimuthal  point.  If  this  point 
subtends  an  angle  8  to  the  array  axis  then  from  Fig. 1  the  electrical 
signal  received  at  subarray  one  is  advanced  by  a  phase  0  where 

0  =  (kd  sin  B)/2  (1) 

>  where  i-.  is  the  wave  number  for  the  transmitted  frequency. 

The  signal  received  at  subarray  two  is  retarded  by  0  .  This 
is  illustrated  in  Fig.  2.  Thus  in  the  absence  of  noise  the  measured 
phase  difference  0j  is  given  by 

0d  =  20 

=  kd  sin  B  (2) 

As  the  signal  received  from  the  array  is  contaminated  by 
noise  the  measured  signals  received  at  the  array  are  given  by 


= 

Si 

+ 

nl 

and 

11^2 

s 

s2 

+ 

n2 

where  mj  and  m2  have  phase  angles  0j  and  02  and  nj  and  n2 
are  the  effective  noise  vectors  at  each  array.  The  angle  between 
nj  and  Sj  is  assumed  to  be  ^  and  n2  and  S2  to  be  t(i2  . 

Now 

tan  (0i  -  6)  ■  n1  sin  (S^  +  n.^  cos  ij^)  (3) 


4 

UNCLASSIFIED 


Fig.  2.  The  relationship  of  m  to  n  and  s. 
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and  =  (Si  +  n±  cos  i{»  i)2  +  (ni  sin  ij<  t)2  (4) 

for  i  =  1,  2,  si  =  I  S±  I  and  n±  =  | ti  | 

In  n  is  not  too  large  (3)  reduces  to 


A6i  =  01  -  e 


n^  3in(\|>i)/s  where  S  =  Sj  =  S2 


Then  if  n^  is  assumed  to  be  Gaussian  <|;  uniformly  distributed  and 
E  (*)  denotes  the  expected  value 


E(A0i)2 


(ni2/S2)  (1/2tt) 


■2  TT 


sin2U<.  dto. 
i  i 


(n^/lS2) 


(5) 


(1/2  S  ) 
e 


where  Sg  is  the  expected  signal  power  divided  by  the  expected  noise 
power  at  the  subarray.  Also 


E(m12) 


S2  +  2  (cos2  ip^  +  sin2 


r2w 


+  2  Sn  (l/2n) 


cos  ip  dij^ 


S 1  +  n. 


(6) 


As 

E(A0t) 

— 

0 

and 

E(m1) 

= 

S 

the  standard  deviations  of  A0  and  m  are  given  by 


Sd(A0.)  -  1/ (2S  )** 

1  e 

Sd(mi>  *  nt 


and 
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3.  THE  EFFECT  OF  INDEPENDENT  INTERFERING  NOISE 


Now  if  ni  and  n2  may  be  assumed  independent  then  as  the 
measured  phase  error  is  given  by 

A0  =  A0j  —  A02 


the  standard  deviation  of  the  measured  phase  error  is  given  by 

Sd(A0)  =  l/(Se)1'5  (7) 

A  more  useful  form  of  expressing  (7)  is  to  let  Sx  be  the  signal  to 
noise  ratio  required  to  give  a  standard  deviation  of  Sd(A9). 

Then 

Sx  =  -  20  log  (Sd(A0))  (8) 

Also  the  standard  deviation  in  the  amplitude  found  by  summing 
the  two  sub-arrays  is  given  by 

Sd(mi  +  m2)  =  Jl  n 

As  E(mj  +  m2)  =  2  S 

The  signal  to  noise  ratio  for  the  amplitude  is  given  by 
Sa  =  4S2/  Sd2 (mj  +  m2) 

-  2  S  (9) 

e 

This  represents  3  dB  improvement  in  the  signal  to  noise  ratio  relative 
to  the  signal  to  noise  ratio  at  the  subarray. 

4.  THE  EFFECT  OF  CORRELATED  INTERFERING  NOISE  ON  THE  PHASE  AND 

AMPLITUDE 

Now  if  the  interfering  noise  is  due  to  sea  noise  or  reverberation 
it  may  no  longer  be  valid  to  consider  fij  and  n2  to  be  independent.  To 
analyse  this  situation  it  is  necessary  to  consider  the  distributed  nature 
of  the  noise.  As  illustrated  in  Fig.  3  suppose  that  the  noise  power  at 
angle  8  in  an  element  of  width  Af3  is  given  by  A2(g)AB.  The  noise  field 
may  then  be  represented  by  a  number  of  sources  of  strength 


I 
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Fig.  3.  The  noise  distribution. 
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Ph 

= 

Ah  exp(jYh)AB^  for  h  =  -n/2AB  to  tt/2AB 

(10) 

where 

= 

A(hA6) 

Yh 

= 

y(hA3) 

and 

Then  if 

received 

y(B)  is  the  phase  of  the  noise  generated  at  an  angle  B, 

R(B)  is  the  subarray  response  to  a  noise  source  at  B  the 

from  the  h  th  element  at  subarray  one  has  an  amplitude 

noise 

nih 

= 

Ah  Rh  (A6)^ 

(ID 

where 

«h 

* 

R(hA0) 

and  phase  d  j 

given  by 

dlh 

= 

Y  +  kd  (sin  (hA0)  }  /2 

h 

(12) 

Similarly  at  subarray  two 

Yu  -  kd(  sin  (hAB)  }  /2  (13) 

h 

dlh  -  6 

d2,  +  0 

h 

the  component  in  quadrature  to  is  given  by 

n.,  sin  <l».  for  i  *  1,  2 
in  in 

and  the  component  in  the  direction  of  is  given  by 

n.,  cos  <|>  for  h  *  1 ,  2 
ih  ih 

Thus  if  Sd(A0  )  is  the  standard  deviation  of  the  phase  error  due  to 
n 


T2, 


and 


d  2, 


Then  if. 


+  lv 


and 
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the  hth  noise  source 


Sd(A0  ) 2 
n 

=  E 

{ (A^2/S2)R^2)  {sin  (Yh 

+  kd(sin  6U) /2  -  6) 
h 

- 

sin  ( y  -  kd  (sin  6.  )/2 
n  h 

+  9)}2  AS 

where 

II 

JZ 

CQ 

hAS 

i.e. 

Sd(A0,  )  2 
n 

=  E 

{  (A^2  R^2)/S2}  {4  cos2 

(Yh)  sin2{kd(sin  B^)/ 

2 

-  0}  AS 

(15) 

For  any 

Bh  ’  Yh 

may  be  assumed  to  be  uniformly  distributed  as 

each  range  cell  covers  a  number  of  wavelengths  at  the  frequency  of 
operation. 


Then  as  E  {cos2(yh>}  =  1/2 

Sd(A6h)2  =  2(Ah2/S2)  I^2  sin2  {  kd(sin  &  )/2  -  0}A6  (16) 


Similarly 


Sd (m  ) 2 
ih 


2\2  Rh2 


cos2 


{kd(sin  6,)/2  -  0}  AS 
n 


(17) 


The  square  of  the  standard  deviation  of  the  sum  of  the  components 
Am^^  and  AS^  is  found  by  summing  the  square  of  the  standard  deviation  of 
each  component . 

Thus 


Sd(A0)2  =  l  Sd  (A0h)2 


The  summation  may  be  replaced  by  integration  as  AS  is  made  arbitrarily 
small  so  that 


Sd(A8) 2 


(2/S2) 


r n/2 

A2(6)R2(B)  sin2  {kd  sin  6/2  -  9)  dg 


— TT  /  2 


It  is  now  assumed  that  the  noise  is  isotropic  so  that 


A2 (6)  -  A2 
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and  is  independent  of  6  and  the  array  is  assumed  unshaded  and  continuous 
so  that 


R(B) 


sin  ((kb  sin  g)/2)/((kb  sin  B)/2) 


(18) 


For  kb  >>  1  sin  8  :  8  except  where  R(8)  :  0. 

Then  letting  u  =  k  d/2  and  v  =  k  b/2 

fir/ 2 


Sd(A6)2  =  2(A/S) 2 


sin2(\>B)  sin2(ug  -  9)/(vB)2  dB  (19) 


-tt/2 
•n/2 

-u/2 

cos2  (ug).  sin2  6  }  /(\>g)2  dB 


2(A/S)2 


{sin2(v8)  sin2(u8)  cos2e  +  sin2(vB) 


From 


and 


Q] 

f  n/2 

sin2  (vg)  sin2  (uB)/B2  dB 

-it/ 2 
f  tt/2 

1  sin2  (vB)  cos2  (uB)/B2  dB 
—tt/2 


(tt/2)  {min(u,  v)  } 


(tt/2)  { 2v  -  min(u, v)}  (20) 


Thus 


Sd(A0)2  =  (A2/S2)  tt  (min(u,v)  cos2  0  +  (2v  -  min(u, v))  sin2  0)/v‘  (21) 


Also  the  total  noise  received  by  the  arrays  is  given  by 

(■tt/2 


A2 


sin2  (vB)/(vB)2  dB  =  A2ti/v 


(22) 


—tt/2 


Therefore, 

Sd(A8) 2 


n2/S2  {cos2  0  +  sin2  6} 
n2/S2 

1/S  for  v  £  u 

e 


r 


i 
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and 


Sd(A6)2 


Similarly 


Sd(m)2 


1/Se  {(u/v)  cos2  0  +  (2  -  u/v)sin2  6}  (23) 

for  v  >  u 


n2  for  v  £  u 

n2  {(2  -  u/v)  cos2  0  +  (u/v)  sin2  0}  (24) 

for  u 


The  phase  error  is  minimised  if  0=0  which  may  be  achieved  at  long 
ranges  by  centering  the  target  in  the  beam. 


Then 

Sd(A0) 2 

i.e.  Sd(A0)2 


(d/b  Se)  for  b  >  d 

(1/Sg)  for  b  £  d 

min  (d,b)/b  Sg  (25) 


Similarly 

as 

the  signal 

to  noise  ratio  is 

given  by 

s 

a 

= 

4  S2/(2  Sd(m)2) 

S 

a 

= 

2  S 

e 

if  d  ^  b 

= 

2  Se/  {2  -  b/d) 

if  d  <  b 

i.e. 

S 

a 

= 

2  S  b/(2  b  -  min(b,d)) 
e 

5. 

THE 

AZIMUTHAL 

ERROR 

(26) 


To  determine  the  standard  deviation  of  the  azimuthal  error 
from  the  standard  deviation  of  the  phase  error  it  may  be  noted  in  Fig. 
1  that  if  a  reflecting  point  at  a  range  r  is  a  distance  h  off  the 
array  axis  then 


h  ■  r  sin  g 


(27) 


Then  from  (2) 
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(h/r)  *  0d/kd 

and  thus  Sd(h)  =  (r/kd)  S^C&o) 

Then  substituting  the  value  for  S,(A0)  from  (25) 

a 

Sd(h)2  =  (r2 /k2 d2 ) (min  (b,d)/b  S£)  (28) 

From  (22)  n2  =  A2it/V 

=  27rA2/kb 

Then  S  =  S2 /n2 

e 

=  S2  kb  /  2u  A2  (29) 

Substituting  (29)  in  (28)  gives 

Sd(h)2  =  ( (r2  /k2)  (2irA2/kS2)  (min  b  ,d)  )  /d2b2 

=  q  min  (b,d)/d2b2  (30) 

where  q  =  r2A22ir/k3S2 

Similarly  from  (26) 

Sa  =  2  Sg  b/(2b  -  min(b.d)) 

=  (S2k/irA2)  (b2/ (2  b  -  min(b,d))') 

=  pb2/(2b-  min(b,d)) 

where  p  =  (S2  k/  tr  A2)  (31) 

6.  THE  OPTIMUM  ARRAY  CONFIGURATION 

Now  if  the  total  width  of  the  array  is  c  the  distance 
between  the  subarray  centres,  d  may  be  varied  by  altering  the  sub¬ 
array  width  b  as  shown  in  Fig.  4. 

Then  as  d  *  c  -  b 

Sd(h) 2  -  q  (min(c  -  b),b)/b2(c  -  b)2 
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If  b  <  c/2  then  Sd(h)  is  a  minimum  when 


3 /3b  (l/b(c  -  b)2)  =  o 

or  b  c/3  (32) 

If  b  >  c/2  then  Sd(h)2  is  a  minimum  when 
3/3b  ( 1/b2 (c  -  b))  =  0 

or  b  =  2c/3  (33) 

At  both  these  minima 


Sd(h)2 

= 

6.75  q/c3 

(34) 

From  (31)  if  b 

= 

2c/3 

(35) 

S 

a 

= 

4  p  c/9 

If  b 

= 

c/3 

S 

A 

= 

p  c/3 

(36) 

Thus  of  the  two  solutions  (32)  and  (33)  which  maximise  the  azimuthal 
resolution  setting  b  =  2c/3  gives  a  signal  to  noise  ratio  which  is 
1.2  dB  greater  than  setting  b  =  c/3 

Alternatively  if  optimisation  is  carried  out  for  the  signal  to  noise 
ratio  of  the  amplitude,  then  as 

S  =  p  b2/(2b  -  min  (b,d)) 

For  b  <  d 

S 

a 


p  b2/(2b  -  b) 
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Then  3S  /ab  =  p  (37) 

a 

Thus  the  signal  to  noise  ratio  increases  with  b  for  b  <-  c/2 

If  b  >  d 

Sq  =  p  b2/(2b  -  (c  -  b)) 

=  p  b2/(3b  -  c) 

Then  3S  /ab 

a 

=  p  {  (2b (3b  -  c)  -  b23)/(3b  -  c)2  } 

=  p  {  (3b2  -  2  be) / (3b  -  c)2  } 

=  0  when  b  =  2c/3  (38) 


This  is  the  only  turning  point  in  the  region  c/2  <  b  <  c  and 
corresponds  to  a  minimum.  Thus  Sg  is  a  maximum  in  this  interval  if 


b 

c/2 

or 

b 

c 

(39) 

When 

b 

c/2 

S 

a 

p  c/2 

When  b 

*  c  S  = 
a 

p  c/2 

(40) 

Thus  a  dual  optimum  exists  for  Sg  .  For  b  *  c  however  the  azimuthal 
error  tends  to  infinity.  Thus  if  it  is  required  to  optimise  the  signal 
to  noise  ratio  of  the  array  b  should  be  set  equal  to  c/2. 

Then  Sd  (h)2  -  8  q/c3  (41) 

Relative  to  the  solution  of  (34)  this  gives  a  0.5  dB  improvement  in  the 
signal  to  noise  ratio  for  the  amplitude,  at  the  expense  of  a  degradation 
in  the  phase  performance  which  corresponds  to  a  signal  to  noise  loss  of 
0.7  dB.  The  variations  in  and  Sd(h)  as  b  is  changed  are  illustrated 


.o 


AZIMUTHAL  ERROR  S*(hr  SIGNAL  TO  NOISE  RATIO  (S*) 
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Fig.  5.  The  variation  of  the  phase  error  and  the  signal 
to  noise  ratio  with  Subarray  Width. 
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7.  CONCLUSION 

A  split  beam  sonar  using  unshaded  subarrays  has  been  analysed. 
For  such  a  sonar  the  maximum  azimuthal  resolution  is  obtained  when  the 
two  subarrays  are  each  two  thirds  of  the  width  of  the  full  array  width 
and  the  maximum  signal  to  noise  ratio  is  obtained  when  the  subarrays 
are  one  half  of  the  full  array  width. 
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